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The Hubbard model on a cube was revisited and extended by both nearest-neighbor (nn) 
Coulomb correlation W and nearest- neighbor Heisenberg exchange J. The complete eigensystem 
was computed exactly for all electron occupancies and all model parameters ranging from minus 
infinity to plus infinity. For two electrons on the cluster the eigensystem is given in analytical 
form. For six electrons and infinite on-site correlation U we determinded the groundstate and 
the groundstate energy of the pure Hubbard model analytically. For fixed electron numbers we 
found a multitude of ground state level crossings, in dependence on the various model param- 
eters. Furthermore the groundstates of the pure Hubbard model in dependence on a magnetic 
field h coupled to the spins are shown for the complete U — h plane. The critical magnetic field, 
where the zero spin groundstate breaks down is given for four and six electrons. Suprisingly we 
found parameter regions, where the ground state spin does not depend monotonously on J in 
the extended model. For the cubic cluster gas, i.e. an ensemble of clusters coupled to an electron 
bath, we calculated the density n(/i, T, h) and the thermodynamical density of states from the 
grand potential. The ground states and the various spin-spin correlation functions are studied 
for both attractive and repulsive values of the three interaction constants. We determined the 
various anomalous degeneration lines, where n(/x, T = 0, h = 0) steps more than one, since in 
this parameter regions exotic phenomena as phase separation are to expect in extended models. 
For the cases where these lines end in triple points, i.e. groundstates of three different occupa- 
tion numbers are degenerated, we give the related parameter values. Regarding the influence of 
the nn-exchange and the nn-Coulomb correlation onto the anomalous degeneration we find both 
lifting and inducing of degeneracies in dependence on the parameter values. 



1 Introduction 

The conspicuous renaissance in the study of small Hubbard clusters is mainly due to two rea- 
sons. The first is the availability and success of cluster methods, which were developed within 
the context of strong electron correlation during the last decade [1-6] and applied to problems of 
high-T c superconductivity. A considerable amount of insight into the hot topics of pairing mech- 
anisms, spin-charge separation, charge ordering and pseudo gap behavior has been reached by 
detailed studies of small Hubbard clusters [7-14]. The other source is the technical possibility to 
produce and reproduce high quality nanostructures which makes it possible to study the behavior 
of quantum dot clusters [15-21] or even clusters of a few atoms contingently coupled to organic 
molecules [22-24]. Of course, a detailed knowledge of the cluster physics is inevitable in both 
fields. If one aims at clusters with strongly correlated electrons, the Hubbard model, occasionally 
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extended by nn-Coulomb- and/or exchange interaction, is often a reasonable starting point. Com- 
plete analytical solutions (in the sense, that all eigenvalues and all eigenvectors are determined) 
exist for the Hubbard model with more or less additional terms for the dumpbell, the triangle, 
the square, the tetrahedron, the square pyramid and the octahedron [7,25-30]. Incomplete (but 
very valueable) solutions, e.g. restricted to the ground state or to half or low filling, big external 
fields etc. exist for n-site rings [31,32]. For other models analytical expressions are available for 
special electron (or hole) occupations and/or parameter sets, e.g. for the square with next nearest 
neighbor hopping t' equal to nearest neighbor hopping t. 

The present paper is very related to our former studies of the (extended) Hubbard model on 
three and four-site clusters and the related cluster gases [7,29,30]. It is rather astonishing to 
what a degree it is possible to understand qualitatively the rather involved phase diagram of the 
T c superconductors on the bases of the square-cluster gas model. The main reason for that is 
the degeneration of cluster states with different particle numbers and spin, which has a great 
influence of the coupling of spin- and charge degrees of freedom, and indeed, the degeneracy of the 
ground states for the four-site cluster occupied either with two or four electrons, is responsible for 
the charge separation and inhomogeneities in the pseudogap phase of the cuprates, whereas their 
bosonic character accounts for the superconducting phase [8, 10]. Also the d-wave symmetry of the 
gap seems to be a direct consequence of the symmetry of the four-site cluster ground states [9,13]. 
As was shown in [29,30], the addition of Coulomb or Heisenberg interaction to the pure Hubbard 
Hamiltonian influences this degeneration considerably. 

In the present paper we re-examine the Hubbard model on a cube to study the influence 
of additional nearest-neighbor (nn) Coulomb and nn-exchange interaction, which are added to 
the standard Hubbard model. The main focus of the present paper lies on the grand canonical 
spectrum and its degeneracies between states differing by more than one electron. We want to see, 
where such regions in the complete parameter space of the pure Hubbard model are and how they 
are influenced by the additional interactions. After introducing the model in the next section, 
the third section will be based on a numerical calculation of the complete canonical potential for 
every occupation number. The fourth section is devoted to ensembles of clusters coupled to an 
electron reservoir, the cubic cluster gas. Here, the focus lies on ground state phase diagrams, 
which show the degenerations differing by more than one electron. The behavior in the vicinity 
of the degeneration points or lines may give rise to a likewise rich physics in the "cluster gas" 
or extended systems built from these clusters as it was observed in the relation of the four-site 
Hubbard model on a square to the square lattice. In a concluding section we will discuss the 
results. 

2 The Hamiltonian 

In the following, we consider the model 



H = Hh + He + H j 



(i) 



with H# being the pure Hubbard model: 



u H = t 



(2) 



He is the nn-Coulomb repulsion 




(3) 
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Irreducible representation: I\ T2 T3 T4 T 5 Tq T? Y% Tg Tiq 
Dimension: 1123311233 

Table 1 Dimensions of the irreducible representations of the cubic point group Oh- 
Furthermore a nn-Heisenberg exchange term Hj 

H ^ = \ E s ^ ( 4 ) 

was added. Here c+ and c ia are the creation and destruction operators of electrons at site i with 
spin cr, n icr = c+ c icr and = ^ a n^. indicates the local spin operator at site i. For a detailed 
physical reasoning of the additional terms in the Hamiltonian see e.g. [33]. Our Hamiltonian is an 
extension of the Hamiltonian used e.g. in [34], since we treat the nearest-neighbor correlation W 
and the nearest-neighbor exchange constant J as independent parameters. This is done mainly 
to cover pure Heisenberg type models by setting W/t = 0, J/t ^ 0. The chemical potential [i 
and the magnetic field h in z-direction are introduced to take into account the effects of doping 
and applying external magnetic fields. Please note, that the signs in front of the hopping and 
exchange term are positive, thus one has to be carefully while comparing with papers where 
other conventions are used. There are 4 8 states within that model. Utilizing all symmetries, we 
reduced the Hamilton matrix to a block diagonal form. It has to be noted, that the particle-hole- 
symmetry is broken in case of W/t ^ 0. For details of this symmetry reduction we refer to [7]. In 
contrast to our previous works, the biggest block has dimension 88, what prevented a closed form 
solution. Nevertheless, the smallness of the analytically given block matrices allows a very quick 
and extremely precise evaluation of the eigensystem. The results of the pure Hubbard model were 
compared to the results of Callaway et al. [28] and to a complete independent pure numerical 
calculation of the spectrum [35]. 

One interesting property of the cluster is its equivalence to the simple cubic lattice, using 
periodic boundary conditions in every other bond with kinetic energy set tot/2 [36]. All following 
results are applicable to that system as well. 

3 Fixed electron occupations 

3.1 Rigorous results 

For two electrons on the cluster the Hamilton matrix can be brought to block diagonal form by 
applying all the U-independent symmetries, with the maximum block dimension of four. Thus 
we have for that case both the eigenvalues and eigenvectors as closed analytical expressions. 
We mention here, that an alternate exact solution to the two-electron problem in the domain 
s z — in d-dimensional hypercubes was given before in Ref. [37]. Although we expect that our 
complete solution for the nondegenerate part of the spectrum coincides with their solution in 
the threedimensional case in the considered domain, a comparison was prevented by the very 
different conventions characterizing the eigenstates and a different set of symmetries employed. 
In the appendix we give the symmetries of the groundstates and the corresponding energies in 
dependence on the model parameters in a closed form. 

For the pure Hubbard model (i.e. J=0 and W=0) in the infinite U limit we have investigated the 
Nagaoka state [38,39], which is quite trivial since it has spin | and Ti symmetry. It is an eigenstate 
of the Hamiltonian for finite U also, despite the fact, that it does not depend on the on-site 
correlation. The related energy is — 3t which becomes the groundstate for U > 39.641741191(1). 
For TV = 6 in the infinite U limit the groundstate is highly non-trivial. It has again spin zero and 
Ti symmetry. In contrast to the Nagaoka state it depends monotonously on the on-site correlation. 
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It becomes the groundstate for U > 61.312646262(l)t, where the related energy is — 4.90564t and 
goes to — 2y 4 + y/3t if the on-site correlation goes to infinity. The letter value agrees to 6 digits 
with the numerical calculated number given by Takahashi [40] . 

3.2 Pure Hubbard model 

The pure Hubbard model for the isolated cubic cluster has been discussed convincingly by Callaway 
et al. [28]. We were able to reproduce all figures in the reference and could approve all properties 
mentioned in the text. The ground state spin, given in Table 1 of the reference has been extended 
to negative correlation parameters U. The only difference occurs for n = 5, where the ground 
state has minimal spin S = \ for U < Ot. For all other occupation numbers, the spin for negative 
U is the same as for vanishing U. 

It is due to the rapidly increasing computer power over more than twenty years, that we are 
able now to draw a more detailed picture of the magnetic properties of the system with respect 
to an applied magnetic field h. For better comparison, the complete parameter space is shown in 
the same manner as in [29] using the following mappings: 

u ' := GtTW\ h ' := iTW\ J ' := TTJT\ W ' := TTW\ (5) 

Note, that the correlation parameter has been scaled differently compared to the other parameters, 
mainly to emphasize the transition region between weak and strong correlation, which is located 
around U = 6t, the band width of the simple cubic lattice. 

Although states of high energy might be interesting for e.g. spectroscopy, the present paper 
focuses on the ground states and their dependence on the various parameters. The quantum 
numbers of the lowest states for each set of parameters will be shown in a two-dimensional plot, 
allowing for studying the dependence on two parameters simultaneously. In analogy to thermo- 
dynamics of extended systems, these diagrams are called ground state phase diagrams (GSPDs), 
although they display no phase transitions, in the sense of a diverging correlation length. The 
GSPDs utilize colored areas to depict regions of constant ground state quantum numbers. The 
degeneracy of the respective ground states may be calculated by multiplying the spin degeneration 
25+1 by the spatial one, which is equal to the dimension of the irreducible representations printed 
in Table 1. If a magnetic field is applied, the spin degeneracy is lifted of course. 

The first GSPDs in Fig. 1 analyze the isolated cubic cluster with respect to the on-site correla- 
tion U and an external magnetic field h for all occupation numbers. They obviously are symmetric 
with respect to the magnetic field, as expected. Additionally, the case g) n = 7 is identical to 
i) n = 9 as is required by the particle-hole-symmetry. Since coding errors usually destroy this 
symmetry, we calculated the pictures independently, thus increasing the confidence in our results. 
Nevertheless, in what follows only the bottom half of the occupation numbers will be shown, if 
particle- hole-symmetry holds, what is the case for W/t = 0. 

Even occupation numbers exhibit a non-magnetic ground state for absent magnetic field, which 
is in accordance to the picture of an antiferromagnetic state. Though, it is an interesting feature, 
that in some cases, e.g. n = 4 and n = 6, a very small magnetic field is sufficient at small on-site 
correlation to produce a magnetic ground state, whereas in the other cases n = 2 and n = 8 
no such behavior may be seen. For the interesting cases n = 4 and n = 6, Fig. 2 shows the 
critical magnetic field he destroying the non-magnetic ground state in dependence on the on-site 
correlation, e.g. the line between the white and light red region of figures d) and f). The states 
with S z = are most stable in an intermediate region of U. Similar behavior has been found 
for the tetrahedron at n = 4, but not for the hexagon, showing that this is not generic for small 
clusters. The kink in the n = 6 image is caused by the symmetry transition from a two fold 
degenerated T3 groundstate with spin zero below U = 61.313t to the singlet Ti above [28]. In 
case of an odd occupation number, the picture is much more uniform, in that there cannot be any 
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Fig. 1 Total spin projection S z of the canonical ground states of the isolated cubic cluster in dependence 
on the on-site correlation U and the magnetic field ft for all occupation numbers in increasing order, from 
a) n = 1 to i) n = 9. The colored areas denote a constant value of S z with the inscribed value. To show 
the complete parameter range, the axes have been transformed non- linearly using eq. (5). 
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Fig. 2 Critical magnetic field 
he, where the ground states with 
S z = and —1 are degener- 
ated, in dependence on the on-site 
correlation U. The kink on the 
right panel is caused by the spa- 
tial transition at U — 61.313£. 



states with vanishing S z and the magnetization rises with either increasing magnetic field strength 
| ft | or the on-site correlation U. At fixed U, the system runs through all possible spin projection 
values S z with increasing strength of the magnetic field until saturation. 
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Fig. 3 Canonical GSPDs in dependence on 
U and spin exchange parameter J for the oc- 
cupation numbers n — 2 ... 8. Colored areas 
stand for ground states of constant quantum 
numbers, which are listed on the right of each 
figure. Because no magnetic field is applied, 
all states are degenerated with respect to S z - 
Degenerations may be read off by multiplying 
the spin degeneracy by the dimension of the 
irreducible representation IV 
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Fig. 4 Canonical GSPDs in dependence on 
the on-site (U) and nearest-neighbor coulomb 
interaction (W) for the occupation numbers 
n — 2 ... 8. Colored areas stand for ground 
states of constant quantum numbers, which 
are listed on the right of each figure. Because 
no magnetic field is applied, all states are de- 
generated with respect to S z . Degenerations 
may be read off by multiplying the spin de- 
generacy by the dimension of the irreducible 
representation IV 
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Fig. 5 Canonical GSPDs for the parameters 
of the extended model at U = 6t for occupa- 
tion numbers n = 2 ... 8. The dependence 
on the nearest-neighbor coulomb interaction 
W and the exchange interaction J is shown. 
Colored areas stand for ground states of con- 
stant quantum numbers, which are listed on 
the right of each figure. Because no magnetic 
field is applied, all states are degenerated with 
respect to S z . Degenerations may be read off 
by multiplying the spin degeneracy by the di- 
mension of the irreducible representation IV 
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Fig. 6 Ground state spin S in dependence on the exchange parameter J for the canonical case n — 6 at 
three different on-site correlations U/t — 2,4,8. Note the intermediate region around U = 4t, where the 
spin does not increase monotonically with lower J, what could have been expected and is retrieved for 
the other cases. 



3.3 Extended Hubbard model 

The additional parameters of the extended model introduce a wide variety of new level crossings, 
which might be understood by the help of GSPDs. In Fig. 3 the dependence on the on-site 
correlation U and the spin exchange interaction J is shown for occupation numbers n = 2 to n = 8. 
The cases n < 2 are trivial and therefore omitted. According to the Hamiltonian (4), one has to 
distinguish between a ferromagnetic (J/t < 0) and an antiferromagnetic (J/t > 0) interaction. 
In the first case, the additional term favors spin alignment. For extremely strong ferromagnetic 
interaction, the total spin is increased to its maximum for all U at every occupation number. In 
some cases, e.g. n = 4 or n = 6 even a very small J/t < introduces a change to S = 1 for small 
positive on-site correlations, whereas for n = 5 a small positive J lowers the spin from S = § to ^. 
For the three occupation numbers n = 4, 5, 6 the exchange parameter therefore has a big influence 
on the ground state spin. Contrary to the case n = 5, an antiferromagnetic interaction has no 
drastic effects in general, since the correlation parameter U favors antiferromagnetic alignment, 
too. This is especially visible in the case of the repulsive model, where only a few spatial transitions 
are visible. For n = 7 the afore mentioned Nagaoka state [38,39] is obtained for U > 39.642t for 
the pure Hubbard model. The n=7 panel of Fig. 3 shows the influence of the additional exchange 
interaction on the Nagaoka state (blue): Whereas the maximum spin state is stabilized for J/t < 0, 
it is completely destroyed for approximately J > 0.06t even with U — > oo. 

This influence of the exchange term holds in general, although there are some regimes, where the 
effects are more subtile. Fig. 6 shows the ground state spin S for different on-site correlations over 
a relatively small region from J = — t to J = 0£, where for lower values an increasing tendency for 
ferromagnetic alignment is expected. Indeed, this is found for both weakly and strongly correlated 
systems at U = 2t and U = 8t, respectively. Contrary, for an intermediate region around U = 4t 
the spin does not depend monotone on J showing the complex nature of the interaction. 

The GSPDs for the nearest-neighbor coulomb interaction W are shown in Fig. 4. Interestingly 
enough, there are differences in the effects of both coulomb interactions, most visible in the right 
column, where the total spin increases with high positive values of U but not with W. The 
situation is different for negative U, where no change in total spin is observable, nevertheless 
spatial transitions take place. For the even occupation numbers shown in the left column, the 
total spin is mostly zero. The only exceptions are for n = 4 a region for small W and positive 
U and for n = 6 at small positive U and negative W. Consequently, adding a small nearest- 
neighbor Coulomb interaction rarely introduces ground state changes, notably in the stability of 
the Nagaoka state at n = 7. The picture of the extended model is completed by the GSPDs in 
dependence on the two extended parameters W and J for fixed correlation strength U = 4t shown 
in Fig. 5. The images show the rich interplay of both interactions. Although it is hard to detect 
general rules, one may roughly say that nearest-neighbor Coulomb interaction acts likewise as 



12 



Schumann and Zwicker: Extended Hubbard model on a cubic clusters 




Fig. 7 The occupation number in dependence of the chemical potential \i in the cubic cluster gas for 
three different correlation strengths JJ It — 1, 4, 8 (f.l.t.r.). The two pictures on the left visualize the weakly 
correlated region, whereas the right one shows a strongly correlated system. The dotted curves for a finite 
temperature show the usual smearing of the steps. 




Fig. 8 The occupation number of the ground state in dependence of the chemical potential \i and the 
correlation parameter U in the cubic cluster gas. The left panel shows the complete parameter space using 
primed values, whereas the right one is a magnification with unsealed parameters, to show the small area, 
where the black marked region n — 7 is the ground state. 



the on-site correlation U does. The latter is deduced from the comparison to Fig. 3, where the 
pictures show similar structure. 

4 Cluster gas 

4.1 Cluster gas approximation 

In contrast to the last chapter, where the results for fixed occupation number, were given, the 
present one is devoted to the so called "cluster gas", what is an ensemble of non-interacting 
identical clusters coupled to a reservoir of electrons, thus introducing the chemical potential \i. It 
serves as an approximation to the full extended lattice, where every second bond is replaced by the 
indirect hopping via the particle bath, thus allowing for fluctuations of the occupation numbers 
on each cluster and therefore modeling the effect of doping. 

4.2 Pure Hubbard model 

In a first step, the influence of the chemical potential on the pure Hubbard model will be analyzed. 
The occupation number n in dependence on \i for different values of the on-site correlation is 
shown in Fig. 7. The step functions show similar behavior as in other small clusters [29], e.g. 
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Occupation n 


Correlation UTp/t 


Chem. potential /J^Tp/t 


2,3,4 


13.29885 


-0.4993926 


6,7,8 


3.372781 


0.5218219 


6,7,8 


7.106821 


1.3549741 


6,7,8 


15.32480 


1.8342337 



Table 2 Ground state triple points 
of the cubic cluster gas for U/t > 
and n < 8. The related parameters 
have been calculated to the printed 
accuracy. 




Fig. 9 Lowest energy 
levels in the vicinity of a 
degeneration point in de- 
pendence on ji at U = At. 
A magnetic field h/t = 
0.5 has been added in the 
right panel. 



there are steps higher than one. These steps correspond to degeneration points, where ground 
states differing in their electron occupation by more than one have the same energy. 

Fig. 8 shows an overview of this essential feature. It generalizes the results to cover the complete 
parameter space of both the correlation parameter and the chemical potential in a GSPD using the 
scaling functions (5). Comparing the two figures, it is convincing that this non-linear scaling does 
not change the qualitative picture of the phase diagram and may be used in the following. The 
above mentioned steps reoccur as the boundaries between the colored areas in both images. The 
steps higher than one may be found at borders of two regions, which are not consecutive in the 
occupation numbers displayed in the legend to the right. Additionally, there exist points, where 
degeneration lines intersect, e.g. three different ground states coexist. Those ground state triple 
points may be calculated with high accuracy. The results for the repulsive model are printed in 
Table 2 and may be easily checked versus Fig. 8. They may mainly serve as a reference calculation 
for other theoretical analysis, since it is very unlikely to have those special conditions in any real 
experiment. In an week interacting extended system we have to expect the interplay of three 
quantum phase transitions in the vicinity of the triple points. Regarding the degeneration lines, 
we see, that they are present through a wide parameter range and are therefore reachable by 
experiments, e.g. using electron or hole doping or pressure. 

Since the degeneration points - and therefore especially the highly degenerated ones - are 
considered to be the key to the complexity of the phase diagram, the stability with respect to 
the remaining parameters is of major interest. Fig. 9 shows the destruction of such a point by 
applying a small magnetic field. The states with even occupation are not magnetic, while the 
state with n = 3 is a doublet analogous to the canonical cases in Fig. 1. On the other hand, new 
degeneration points may be introduced by applying an external magnetic field, which can be seen 
in Fig. 9, too, where the n = 2,4 degeneracy is lifted by the field. 

The influence of the magnetic field is studied comprehensively by means of the GSPDs shown 
in Fig. 10, which show a multitude of new ground state degeneration lines and triple points. 
In contrast to the canonical cases, the system does not pass all possible spin projections with 
increasing magnetic field, but rather exhibits large jumps, e.g. at ji' « — 0.6t from S z = to 
\S Z \ = 2. This is accompanied with a change in the occupation number from (n) = 2 to (n) = 4 
and in the spatial symmetry from Ti to Tq, what may be deduced from the figures on the left. The 
system switches abruptly from two electrons with antiparallel spins to a configuration, where four 
spins are aligned. This example shows the rich interplay of charge and spin degrees of freedom 
under the influence of a magnetic field, which is a completely new feature of the cluster gas and was 
not present in the canonical case in Fig. 1. At small magnetic fields, the right figure also displays 
an antiferromagnetic ground state, with mostly minimal spin projection. The only exception is 
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Fig. 10 Occupation number n, spin projection 
S z and spatial symmetry Y of the ground states 
in dependence on the chemical potential fi and 
the magnetic field h at fixed correlation strength 
U — 4£. The meaning of the colors are visu- 
alized in the legend right to each figure. Since 
all pictures cover the same parameter space, the 
areas are related, such that n, S z and T may 
be extracted for any given point. The size of the 
different areas may not be compared directly, be- 
cause of the non-linear transformed parameters. 



the case (n) = 5 around h = \i = Ot, where four spins seem to be aligned, while the remaining one 
is antiparallel. This is in accordance with the related canonical results. 

The spatial order of the local spins can be discussed in great detail using spin correlation 
functions, what was first done for the canonical case by Callaway et al. [28]. The grand canonical 
case is shown in Fig. 11. The degeneration lines of Fig. 8 are once again visible, since the same 
parameter range is covered and the spin configuration depends heavily on the number of electrons. 
Consequently, the mean occupation numbers shown in Fig. 8 may be used to explain features of 
the correlation functions. 

The key feature is visible in the top right picture, where the nearest-neighbor correlation is 
shown: Only the region with (n) = 7 exhibits a positive spin correlation denoting a tendency 
for a ferromagnetic ground state, which is in accordance with Nagaoka's theorem [38,39], since 
(n) = 7 is the occupation number, where exactly one electron is missing to half filling. Moreover, 
the calculation shows, that the Nagaoka state is not only stable for infinite on-site correlation 
U, but for the whole region (n) = 7 with U > 39.642t, where the transition to the S = | state 
takes place. An interesting conjecture may be drawn regarding the spatial change of the state 
(n) = 6 at U = 61.313t, which occurs in the vicinity of the Nagaoka state. The switch of the 
irreducible representation from T$ to Ti may be induced by the ferromagnetic n = 7 state, which 
has Ti-symmetry, too. 

All other regions have a tendency for antiparallel spins, leading to antiferromagnetism in the 
extended lattice. The next pictures support this point of view, since especially for half filling 
(n) = 8 parallel alignment is favored across the face diagonal in c), but once again reversed across 
the space diagonal in d), which is in accordance with the view of an alternating spin configuration. 
Contrary, in the case of (n) = 7, all correlation functions are positive, supporting a parallel order. 
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Fig. 11 The four different spin correlation functions of the cubic cluster gas in dependence on the 
correlation parameter U and the chemical potential fi. The values of the functions are visualized by 
colors, which might be associated using the legend right to each figure. The four functions shown are the 
same as in Ref. [28]: Picture a) shows the mean local spin square (Sz,i), whereas the other cases correlate 
spins on different sites i and j, with increasing distance to each other: (S z ,iS z ,j)> There are three distinct 
possibilities on the cube: b) nearest-neighbor, c) face diagonal and d) space diagonal. 




Fig. 12 Thermodynamical density of states for the cubic cluster gas in dependence on the chemical 
potential for different correlation parameters at the temperature k&T — O.lt. The function is shown for 
three increasingly strong correlations U/t = 0, 4, 16. The axis of the chemical potential \x has been shifted 
by U/2t to emphasize the particle- hole-symmetry. 



The picture is less clear for the smaller occupation numbers, since the spatial configuration of the 
spins seems to be more complex and a connection between neighboring spins can not be drawn 
that easily. This is also supported by the overall smaller absolute value of the correlation functions, 
showing the weaker correlation directly. 

To complete the picture of the cubic cluster gas, this section ends with some further character- 
istic functions, which were also computed for e.g. the square [411. The spectral function J + 
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of the cubic cluster has already been reported in Ref. [42] and has been confirmed. Additionally, 
the thermodynamic density of states shown in Fig. 12 is defined as 

With increasing on-site correlation the initial four-peak structure of the uncorrelated model trans- 
forms into a more complex picture: Two groups of peaks separate with an increasing gap inbe- 
tween. The separation takes place at the transition region around U = 6t. This is similar to the 
spectral function and in qualitative accordance with the common picture of the Hubbard-Mott 
transition. 



4.3 Extended Hubbard model 

In the same manner as in the canonical case, additional interactions may be added for the cluster 
gas. To compare the results to the canonical case, the same terms are considered. It is evident, 
that the complete phase diagram is extremely complex and only a small part will be discussed 
here. 
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Fig. 13 GSPDs of the extended model of the cubic cluster gas. The correlation parameter U = At is fixed 
in both figures. The groundstates are pictured as colored areas in dependence on the chemical potential 
fi and both additional parameters J and W, respectively. The quantum numbers are shown only for large 
areas in the legend right to the figures. 



The influence of a ferromagnetic ( J < 0£) and antiferromagnetic (J > Ot) exchange interaction 
is shown on the left panel of Fig. 13. The multitude of degeneration points increases drastically. 
Qualitatively, the same features as in the canonical case of Fig. 3 are visible, in that a ferromagnetic 
exchange leads to charge ordering and the opposing antiferromagnetic one does not have those 
qualitative impact, but stabilizes the tendency for antiparallel spin ordering. Nonetheless, both 
characteristics lead to a half-filled state for extreme values showing the tendency for an ordered 
state in general. The degeneration points discussed in section 4.2 are mostly visible throughout 
the whole range of J, although the state n = 7 is stable around J « —t. New degeneration points, 
e.g. between n = and n = 2 are introduced for positive J. 

The nearest-neighbor Coulomb interaction is shown of the right side of the same figure. The 
qualitative picture is quite similar to Fig. 8, the GSPD in dependence on U and \i. This is 
reasonable considering both U and W denote Coulomb interactions, only differing in their spatial 
properties. More quantitatively, it has been shown for the triangle, the tetrahedron and the square, 
that both parameters occur in fixed combinations in most eigenvalues of the problem [29,30]. The 
basic degeneration points of section 4.2 are stabilized for high absolute values of W and new points 
are introduced, where n steps e.g. from 4 to 8 at W = t. 
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5 Discussion 

The comparison of our semi-analytical treatment of the Hubbard model on a cubic cluster with 
the former results of our analytical solutions on smaller clusters makes clear, that the complexity 
is increased considerably by increasing the cluster size or by adding additional interactions to the 
Hamiltonians. For eight sites we were still able to get the exact spectrum and eigenfunctions for 
the complete parameter region from — oo to +00 for all parameters, thus providing a reference for 
other studies, where the (extended) Hubbard model is employed. 

The detailed reproduction of the results by Callaway et al. is a result, which provides us 
with confidence in our symbolic code. It was a prerequisite for the treatment of the Hubbard 
model extended by additional terms, i.e. an external magnetic field h coupled to the spins, and 
the nearest neighbor Coulomb and exchange interaction, W and J, respectively. The sometimes 
complex effects of increasing h have also been found in other small clusters showing the competition 
between aligning the spins and minimizing their kinetic energy. The situation becomes even more 
difficult with the additional interactions: Although the tendency for ferromagnetic alignment 
caused by J/t < has been recovered for most parameter regions, there are exceptions, where 
lowering J also lowers the total spin, which would have rather been expected for positive J. The 
influence of the nearest neighbor Coulomb interaction W is comparable to the on-site correlation 
U and does mostly introduce spatial transitions in the GSPDs. It is questionable, whether the 
enormous amount of spatial changes may be seen in bulk systems, but it has been recovered that 
the difference between odd and even electron occupation is of qualitative character, which has 
been reported for all other small clusters, too. 

With regard to extended systems, the cubic cluster gas is a model for the simple cubic lattice. 
The approximation will be not too bad, if one assumes a weak inter-cluster hopping, which is 
than replaced by the indirect exchange via the particle reservoir. Besides the highly non-trivial 
analytical closed- form solution for two electrons, the most important results reported in this paper 
are the ground state level crossings being visible in the various GSPDs. They may give rise to 
quantum phase transition in extended systems. With respect to the cluster gas, the degeneration 
points of states with occupation numbers differing by two or more are of central interest, e.g. the 
jump from n = 6 to n = 8 visible in Fig. 7. This is especially important, since this degeneration 
line is present at moderate on-site correlation and lies in the weakly underdoped region, thus being 
reachable by experiments. Utilizing doping, the mean electron occupation of the cluster gas may 
be fixed to (n) = 7 and the result will be a mixture of clusters with six and eight electrons, if U 
lies between 3.37t and 7.10£ or is above 15.32£. On can presume, that something similar to the 
stripe structure of 2D-systems may be observed. 

The alignment of the spins has been studied using the spin correlation functions and it has been 
found, that the half-filled case shows a tendency for antiparallel alignment, as expected. If the 
on-site correlation is high enough, introducing one hole alters the system completely, such that the 
n — 7 spins are parallel aligned, which accounts for a ferromagnetic ground state in the extended 
system, showing that Nagaoka's result is valid down to U = 39.642t. As a consequence one has 
to expect a ferromagnetic groundstate for a mean electron occupation (n) = | and large on-site 
correlation U > 79.284t if the cubic lattice is approximated by an array of week interacting cubic 
clusters with periodic boundary conditions. For n = 6 we have a singlet ground state with Ti 
symmetry for strong electron correlation U > 61.313£ , which may be considered as the two- hole 
counterpart of the Nagaoka state, apart from the fact that the Nagaoka state does not depend on U, 
whereas the n = 6 ground state does. Moreover, the limit value lini£/_ >00 E(n = 6, S = 0, T\)/t = 
-2^4 + V3 may serve as an excellent proving tool for qualifying numerical codes. 

Regarding the influence of the two additional parameters, the nearest-neighbor Coulomb in- 
teraction W and the spin exchange interaction J, it has been found, that the effect of W is 
comparable to the on-site correlation J7, which confirms the finding that the combination of U 
and W may be replaced by an effective on-site correlation parameter in most cases, what was 
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shown analytically for smaller clusters. Contrary, the effect of J is more important, since even 
small values may induce a change of the ground state spin in some areas. An striking example is 
the vanishing of the Nagaoka ground- state in the presence of a small anti- ferromagnetic exchange. 
Additionally, the simple picture of J emphasizing the tendency for (anti-) parallel spin alignment 
does not always hold, since there are parameter regions, where the contrary has been observed. 

The qualitative consequences of the degeneration points onto extended systems was discussed 
in detail in [29] for smaller clusters. Since our results in general fit well to the given scenario, we 
abstain from discussing it again. One remarkable new feature is the existence of two separated 
regions of the correlation parameter U where a 6 — 8 degeneration happens. It is clear that we 
may expect the most interesting effects in the vicinity of these degeneration points, i.e. for n = 7 
or a hole density of 0.125 %. 

Regarding frustration, our results fortify the view, that it is impossible to classify a cluster as 
frustrated from geometry alone. It is inevitable to consider it together with the electron occupation 
and the correlation strength. Typical frustrated situations will cover only small areas in the U-fi- 
GSPDs or even be absent. On the contrary the area of unfrustrated configurations will grow at the 
costs of the area of neighboring frustrated configurations. This reflects the fact, that frustration 
increases the ground state energy, thus only a small amount of chemical energy A/i((7V)) is needed 
to suppress it. 

Probably, it will be hard, to realize an experimental situation, where strongly correlated elec- 
trons reside on a cubic cluster or simple cubic lattice. Most probably seems to be a situation were 
d- or f-atoms are embedded in organic molecules. An other possibility would be the realization 
of a compound, were an effective simple cubic model remains after projecting out other parts of 
the full Hamiltonian, in analogy to the reduction of Hamiltonian of the copper oxides planes to 
an effective Hubbard model on a 2-d simple cubic lattice. 

A final conclusion from our work on small clusters will be in order. We have meanwhile the 
complete eigensystem for all clusters, where the correlation independent symmetries are enough to 
break down the Hamilton matrix to blocks lower than 5. Sorry to say the cube does not belong to 
that set. Nevertheless, the big number of symmetry operators made it possible to block-diagonalize 
the model analytically, with block sizes allowing a diagonalization to arbitrary numerical accuracy 
by help of symbolic computer facilities. This exact solution allowed us to get an impression of 
the unforeseen complexity of the model. We therefore do not intend to go to bigger cluster sizes 
or further additional terms in the Hamiltonian, since with increasing size of the final blocks we 
will loose the benefit from our analytical treatment, therefore, ab initio numerical algorithms will 
be more appropriate for such tasks. We see the main benefit from our analytical/exact cluster 
solutions in its reference function for other methods dealing with strongly correlated electrons and 
as toy models for problems intimately connected to small clusters, e.g. entanglement, quantum 
dot arrays or electron transport through single molecules. 

Acknowledgements We are much indebted to J. Richter for his lively interest in the subject and 
valuable discussions. 
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A Analytical solution for N=2 

In the following we give the eigensystem for the electron occupation N = 2. The original Hamil- 
ton matrix in that space is of dimension 120, which is symmetry reduced to 61 subspaces with 
maximum dimension 4. Thus, it is easy to get the eigensystem in its explicite dependency on the 
correlation parameters Unfortunately the lengthy form of some explicite eigenvalues and eigen- 
vectors prevents printing them. Here we give the closed form expression of the groundstate for 
anti- ferromagnetic or small ferromagnetic exchange (gray area in Figs. 3 and 4) belonging to the 
non-degenerate eigenstate with T\ and 5 = 0: 



Eiq-GS = - 




A 9 Aio 2A 12 ^ 

1 ~ 8 A^M 32 2 /3^ + 3 6^2 

A = A 14 ^10 ^12 ¥M 

2 16 + 32 2 /3^+ 3 6^2 
A 3 = A 7 + V^U 

A, = Al-AAl 

A 5 = -2A\ 2 + 1728A n A 12 - 36A 13 A 14 _A 12 + 864A 2 3 + 324A 11 A 14 

A 6 = A\ 2 + 288An + 12A 13 A 1A 

A 7 = -2A 3 12 + 1728AiiAi 2 - 36A 13 A 14 _A 12 + 864A^ 3 + 324AnA 2 4 

A 8 = A\ 2 + 288An + I2A13A14 

A 9 = ^12^14 + 16^3 

o 

A w = 324JV + 9J 2 U 2 + 624Jt 2 U - 1728Jt 2 W - 48JU 2 W 

+13312t 4 + 336t 2 C/ 2 - W64t 2 UW + 22>0U 2 W 2 + MU 2 W 2 
A n = 3Jt 2 U + 24i 4 - 8t 2 UW 
An = 3JU + 80t 2 - 8UW 
A u = 9Jt 2 - Ut 2 U - 24t 2 W 
A 14 = 3J-2U-8W 

and the groundstate for strong ferromagnetic exchange (dark yellow area in Fig. 3) is ninefold 
degenerated, since it belongs to Tg and s = 1. The related energy is 



1 _i / (J + SW) (J 2 + 16JW - UAt 2 + 6AW 2 ) \ 
B ' " 3 C0S " { ^JUp J 

A l = -L(J + 8H') 2 + 16( 2 . 

Therefore we decided to give the analytical form of the symmetry reduced Hamilton matrices 
together with the basis, utilized for the calculation. This form contains the full information in 
a more condensed form, since the coefficients of the basis are independent on the model param- 
eters up to this stage of diagonalization. The headline to every subspace gives the eigenvalues 
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of the J7, W, J-independent symmetry operators, i.e. the eigenvalues of S together with the spin- 
projection S z in z-direction m s . Furthermore, the spatial symmetry is indicated by Tij, where 
the first index labels the irreducible representation of the point group and the second numbers 
the partner. The notation is based on Ref. [43]. Next we give the basis, where we have chosen the 
same notation as in Ref. [7] and finally the related Hamilton block matrix. 



Subspace No. 1 with N = 2, s = 2, m s = — 1, T 4 _i 



|^l-l) = ~^J2 (l 000000d ^ _ |0000dd00) — |000d000d) + |00d000d0) 
— |00dd0000) + |0d000d00) — |d000d000) + |dd000000)) 
|^l-2) = — 7= (|00000d0d) — |0000d0d0) — |000dd000) + |00d00d00) 

2y 2 

-f |0d0000d0) — |0d0d0000) — |d000000d) + |d0d00000)) 



^ = ' 2 2t 



Subspace No. 2 with N = 2, s = 2, m s = -1, T 4 _ 2 



|^2-l) = -7= (|00000d0d) + |0000d0d0) — |000d00d0) — |00d0000d) 

2y 2 

-f |0d00d000) — |0d0d0000) + |d0000d00) — |d0d00000)) 
|^2-2) = -7= (|00000dd0) + |0000d00d) — |000d000d) — |00d000d0) 

2y 2 

-f |0d000d00) — |0dd00000) + |d000d000) — |d00d0000)) 

_ , 2t 
; ' " 1 2t \(J + %W) 



Subspace No. 3 with N = 2, s = 2, m s = —1, r 4 _ 3 



|^3-l) = —7= (|000000dd) + |00000dd0) — |0000d00d) + |0000dd00) 

2y 2 

-f |00dd0000) + |0dd00000) — |d00d0000) + |dd000000)) 
|^3_2) = \= (|000d00d0) — |000dd000) — |00d0000d) + |00d00d00) 

2y 2 

— |0d0000d0) + |0d00d000) + |d000000d) — |d0000d00)) 

Ha = 1 2t 
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Subspace No. 4 with N = 2, s = 2, ra s = —1, T 5 _i 



|^4_l) = — ^ (|000000dd) — |00000dd0) + |0000d00d) + |0000dd00) 

2y 2 

-f |00dd0000) — |0dd00000) + |d00d0000) + |dd000000)) 
|^4_2) = \= (|000d00d0) + |000dd000) — |00d0000d) — |00d00d00) 

2y 2 

-f |0d0000d0) + |0d00d000) — |d000000d) — |d0000d00)) 

_ , \{J+W) 2t 



Subspace No. 5 with N = 2, s = 2, ra s = —1, T5_2 



|^5-l) = —7= (|000000dd) — |0000dd00) + |000d000d) — |00d000d0) 
2y 2 



— |00dd0000) — |0d000d00) + |d000d000) + |dd000000)) 
— (|00000d0d) — |0000d0d0) + |000dd000) — |00d00d00) 
— |0d0000d0) — |0d0d0000) + |d000000d) + |d0d00000)) 



u _ / W + «W) 2t 
5 1 2t 



Subspace No. 6 with N = 2, s = 2, ra s = —1, r 5 _ 3 



|^6-l) = -7= (|00000d0d) + |0000d0d0) + |000d00d0) + |00d0000d) 

2y 2 

— |0d00d000) — |0d0d0000) — |d0000d00) — |d0d00000)) 
|^6-2) = -7= (|00000dd0) + |0000d00d) + |000d000d) + |00d000d0) 

2y 2 

— |0d000d00) — |0dd00000) — |d000d000) — |d00d0000)) 

2t 



Hn ~ 1 2t ±(J + 8W) 
Subspace No. 7 with TV = 2, s = 2, m s = —1, T 7 



|^7-l) = \= (|000000dd) — |00000dd0) + |0000d00d) + |0000dd00) + |000d000d) — |00d000d0) 

2y3 

— |00dd0000) + |0d000d00) + |0dd00000) — |d000d000) — |d00d0000) — |dd000000)) 
|^7_ 2 ) = — ^ (|000d0d00) — |00d0d000) + |0d00000d) — |d00000d0)) 

l(J + 8W) 
H? ! 
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Subspace No. 8 with N = 2, s = 2, m s = — 1, T 8 _i 



l^s-l) = -7= (|000000dd) — |00000dd0) + |0000d00d) + |0000dd00) 

2y6 



-|00dd0000) + |0dd00000) — |d00d0000) — |eM)00000)) 

z (|000d000d) 
) 

H 8 = ( i(J + 8W) ) 



H — (|000d000d) — |OOdOOOdO) + |OdOOOdOO) — |dOOOdOOO)) 

V6 



Subspace No. 9 with N = 2, s = 2, m s = -1, r 8 _ 2 



|^9-l) = -7= (|000000dd) + |00000dd0) — |OOOOdOOd) + |OOOOddOO) 

2y 2 

— |OOddOOOO) — |OddOOOOO) + |dOOdOOOO) — |ddOOOOOO)) 

H 9 = ( i(J + 8W) ) 



Subspace No. 10 with N = 2, s = 2, m s = —1, T 9 _i 



|^10-l) = — - 7= (|00000d0d) + |0000d0d0) — |000dd000) — |00d00d00) 

2y 2 

-f |0d0000d0) + |0d0d0000) + |d000000d) + |d0d00000)) 
1^10-2) = ^ (|00000dd0) + |0000d00d) + |0dd00000) + |d00d0000)) 

1^10-3) = — ^ (|000d0d00) + |00d0d000) — |0d00000d) — |dOOOOOdO)) 

2y/2t 2y/2t 

H10 = ( 2y/2b \(J + SW) 
2y/2t 



Subspace No. 11 with N = 2, s = 2, ra s = —1, r 9 _i 



|^ll-l) = — ^ (|000000dd) — |0000dd00) + |00dd0000) — |dd000000)) 
I^H-2) = -7= (|00000d0d) — |0000d0d0) — |000d00d0) + |00d0000d) 

2y2 

-f |0d00d000) + |0d0d0000) — |d0000d00) — |d0d00000)) 
|^H_3) = ^ (|000d0d00) — lOOdOdOOO) — |0d00000d) + |d00000d0)) 

±(J + 8W) 2>/2t 
H n = ( 2y/2t 2>/2t 

2V2t 
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Subspace No. 12 with N = 2, s = 2, m s = -1, r 9 _ 3 



|^12-l) = - (|000d000d) + |00d000d0) + |0d000d00) + |d000d000)) 
1^12-2) = — - 7= (|000d00d0) + |000dd000) + |00d0000d) + |00d00d00) 

2y2 

-f |0d0000d0) + |0d00d000) + |d000000d) + |d0000d00)) 
1^12-3) = ^ (|000d0d00) + |00d0d000) + |0d00000d) + |d00000d0)) 

\(J + W) 2V2t 
H i2 = I 2y/2t 2y/2t 



2y/2t 



Subspace No. 13 with N = 2, s = 2, ra s = —1, Tio-i 



|^13-l) = —7= (|000d00d0) — |000dd000) + |00d0000d) — |00d00d00) 

2y2 

— |0d0000d0) + |0d00d000) — |d000000d) + |d0000d00)) 

H13 = ( ) 



Subspace No. 14 with N = 2, s = 2, m. s = —1, Tio-2 



|^14-l) = — 7= (|00000d0d) + |0000d0d0) + |000dd000) + |00d00d00) 

2y 2 

— |0d0000d0) + |0d0d0000) — |d000000d) + |d0d00000)) 

H14 = ( ) 



Subspace No. 15 with N = 2, s = 2, m, s = —1, Tio-3 



|^15-l) = -7= (|00000d0d) — |0000d0d0) + |000d00d0) — |00d0000d) 

2y 2 

— |0d00d000) + |0d0d0000) + |d0000d00) — |d0d00000)) 

H15 = ( ) 
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Subspace No. 16 with N = 2, s = 0, m s = 0, Ti 



|*16-l) 
1*16-2) 



1*16-3) = 



I* 



16-4/ 



Hi6 = 



2^2 



(|00000002) + |00000020) + |00000200) + |00002000) 



1 



2V6 



+|00020000) + |00200000) + |02000000) + |20000000)) 

(|000000du) — |000000ud) + |00000du0) — |00000ud0) + |0000d00u) + |0000du00) 



1 



— |0000u00d) — |OOOOwdOO) + |000d000u) — |000u000d) + |OOdOOOuO) + |OOduOOOO) 

— |OOuOOOdO) — |OOudOOOO) + |OdOOOuOO) + |OduOOOOO) — |OuOOOdOO) — |OudOOOOO) 
-f |dOOOuOOO) + |dOOuOOOO) + |duOOOOOO) — |uOOOdOOO) — |uOOdOOOO) — |udOOOOOO)) 

: (|OOOOOdOu) — |OOOOOnOd) + |OOOOdOuO) — |OOOOuOdO) + |OOOdOOuO) + |OOOduOOO) 



2^6 

— |OOOnOOdO) — |OOOndOOO) + |00d0000u) + |OOdOOuOO) — |00u0000d) — |OOwOOdOO) 
-f |OdOOOOuO) + |OdOOnOOO) + |OdOuOOOO) — |OuOOOOdO) — |OuOOdOOO) — |OnOdOOOO) 
-f |dOOOOOOu) + |dOOOOnOO) + |dOuOOOOO) — |uOOOOOOd) — |uOOOOdOO) — |uOdOOOOO)) 
1 , 

— (jOOOdOitOO) — |OOOnOdOO) + |OOdOuOOO) — |OOuOdOOO) 

2y 2 

-f |OdOOOOOu) — |OnOOOOOd) + |dOOOOOuO) — |uOOOOOdO)) 

/ U -2y/3t \ 

-2V3t AW - ^ At 

At 2y/%b 

\ 2y/Zb / 



Subspace No. 17 with N = 2, s = 0, m s = 0, T 3 _i 



*17-l) = J (|000000dn) — |000000nd) — |OOOOOduO) + |OOOOOudO) — |0000d00u) + |OOOOduOO) 

-f |0000u00d) — |OOOOndOO) + |OOduOOOO) — |OOudOOOO) — |OduOOOOO) + |OudOOOOO) 
— IdOOnOOOO) + |dnOOOOOO) + |nOOdOOOO) — |wdOOOOOO)) 



1*17-2) = 



H 



17 



- (|OOOdOOnO) — |OOOdnOOO) — |OOOuOOdO) + |OOOudOOO) + |00d0000u) — |OOdOOuOO) 

— |00n0000d) + |OOnOOdOO) — |OdOOOOuO) + |OdOOuOOO) + |OuOOOOdO) — |OuOOdOOO) 

— IdOOOOOOn) + |dOOOOnOO) + |nOOOOOOd) — |uOOOOdOO)) 

AW - ^ 2t 



2t 
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Schumann and Zwicker: Extended Hubbard model on a cubic clusters 



Subspace No. 18 with N = 2, s = 0, m s = 0, T 3 _ 2 



|^18-l) = — " 7= (|000000du) — |000000ud) + |00000du0) — |00000ud0) + |0000d00u) + |0000du00) 

— |0000u00d) — |0000ud00) + |00du0000) — |00ud0000) + |0du00000) — |0ud00000) 
-f |d00u0000) + |du000000) — |uOOdOOOO) — |udOOOOOO)) 
1 



r (|000d000u) — |000u000d) + |OOdOOOuO) — |OOuOOOdO) 

2y 3 

-f |OdOOOuOO) — |OuOOOdOO) + |dOOOuOOO) — |uOOOdOOO)) 
1^18-2) = — " 7= (|OOOOOdOu) — |OOOOOuOd) + |OOOOdOuO) — |OOOOuOdO) 

2y 3 

-f |OdOuOOOO) — |OuOdOOOO) + |dOuOOOOO) — |uOdOOOOO)) 
1 , 

(|OOOdOOuO) + |OOOduOOO) — |OOOuOOdO) — |OOOudOOO) + |00d0000u) + |OOdOOuOO) 



4V"3 

— |00n0000d) — |OOnOOdOO) + |OdOOOOuO) + |OdOOuOOO) — |OuOOOOdO) — |OuOOdOOO) 
+|dOOOOOOn) + |dOOOOnOO) — |uOOOOOOd) — |uOOOOdOO)) 



Hl8 = 1 2t 



Subspace No. 19 with N = 2, s = 0, m s = 0, T 5 _i 



|^19-l) = 7= (|00000002) — |00000020) + |00000200) — |00002000) 

2y2 

+100020000) — |00200000) + |02000000) — |20000000)) 

1^19-2) = -7= (|00000d0u) — |OOOOOnOd) — |OOOOdOuO) + |OOOOuOdO) 

2y 2 

-f |OdOuOOOO) — |OnOdOOOO) — |dOuOOOOO) + |uOdOOOOO)) 
|^ig_ 3 ) = (|000d000u) — |000n000d) — |OOdOOOuO) + |OOuOOOdO) 

2y2 

-f |OdOOOuOO) — |OnOOOdOO) — |dOOOnOOO) + |nOOOdOOO)) 
1^19-4) = -7= (|OOOdOuOO) — |OOOnOdOO) — |OOdOuOOO) + |OOuOdOOO) 

2y2 

-f |OdOOOOOu) — |OnOOOOOd) — |dOOOOOuO) + |uOOOOOdO)) 



/ 


u 





-2t 


\ 













2t 




-2t 





AW- ^ 





V 





2t 





/ 
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Subspace No. 

|*20-1> = 
|^20-2) 



20 with N = 2, s = 0, m s = 0, T 5 _ 2 
1 



2a/2 



(|00000002) — |00000020) — |00000200) + |00002000) 



I* 



20-3/ 



2V2 



2V2 



|*20-4> 

H20 = 

Subspace No. 
|^2i-i) = 

|*21-2> 
|*21-3> 
1^21-4) 



2V2 



-|00020000) + |00200000) + |02000000) — |20000000)) 
(|00000dn0) — |00000nd0) — |0000d00u) + |0000u00d) 
-|0dn00000) + |0nd00000) + |d00u0000) — |u00d0000)) 
(|000d00n0) — |000n00d0) — |00d0000n) + |00u0000d) 
-|OdOOnOOO) + |OnOOdOOO) + |dOOOOuOO) — |uOOOOdOO)) 
(|OOOdOnOO) — |OOOnOdOO) — |OOdOuOOO) + |OOuOdOOO) 



— |OdOOOOOn) + |OnOOOOOd) + |dOOOOOuO) — |wOOOOOdO)) 





u 


-2t 





\ 




-2t 





















2t 


V 








2t 


/ 



21 with N = 2, s = 0, m s = 0, r 5 - 3 

~J= (|00000002) + |00000020) — |00000200) — |00002000) 



2a/2 



1 



2a/2 



-|00020000) — |00200000) + |02000000) + |20000000)) 
(|000000dn) — |000000nd) — |OOOOduOO) + |OOOOudOO) 
-|OOdnOOOO) + |OOudOOOO) + |duOOOOOO) — |udOOOOOO)) 
(|OOOdOnOO) — |OOOnOdOO) + |OOdOuOOO) — |00<uOdOOO) 



2a/2 



-|OdOOOOOn) + |OnOOOOOd) — |dOOOOOuO) + |uOOOOOdO)) 
(|OOOduOOO) — |OOOndOOO) + |OOdOOuOO) — |OOuOOdOO) 



|OdOOOOnO) + |OnOOOOdO) — |dOOOOOOu) + |wOOOOOOd)) 



H 



21 





u 


-2t 





\ 




-2t 


4W- ^ 



















2t 










2t 


/ 



Subspace No. 
|^22-i) = 



22 with N = 2, s = 0, m s = 0, T 7 

1 , 

(|OOOOOdOn) — |OOOOOnOd) — |OOOOdOuO) + |OOOOuOdO) — |OOOdOOuO) — |OOOduOOO) 



2VQ 



H 



22 



+ |OOOuOOdO) + |OOOitdOOO) + |00d0000n) + |OOdOOnOO) — |00w0000d) — |OOnOOdOO) 
— |OdOOOOnO) — |OdOOuOOO) — |OdOuOOOO) + |OuOOOOdO) + |OuOOdOOO) + |OuOdOOOO) 
-f |dOOOOOOu) + |dOOOOnOO) + |dOuOOOOO) — |uOOOOOOd) — |uOOOOdOO) — |uOdOOOOO)) 

(0) 
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Schumann and Zwicker: Extended Hubbard model on a cubic clusters 



Subspace No. 23 with N = 2, s = 0, m s = 0, T 8 _i 



|^23-l) = — " 7= (|00000d0n) — |00000u0d) — |0000d0u0) + |0000u0d0) 

2y3 

— |0d0n0000) + |0n0d0000) + |d0u00000) — |u0d00000)) 

H \= (|000d00u0) + |000dn000) — |OOOuOOdO) — |OOOndOOO) — |00d0000u) — |OOdOOuOO) 

4V3 

-f |00u0000d) + |OOnOOdOO) + |OdOOOOuO) + |OdOOuOOO) — |OuOOOOdO) — |OuOOdOOO) 

— |dOOOOOOn) — |dOOOOnOO) + |uOOOOOOd) + |uOOOOdOO)) 

H 23 = ( ) 



Subspace No. 24 with N = 2, s = 0, m s = 0, r 8 _ 2 



|^24-l) = — J (|OOOdOOuO) — |OOOduOOO) — |OOOnOOdO) + |OOOudOOO) — |00d0000u) + |OOdOOuOO) 
-f |00u0000d) — |OOnOOdOO) — |OdOOOOuO) + |OdOOuOOO) + |OuOOOOdO) — |OuOOdOOO) 
-f |dOOOOOOu) — |dOOOOnOO) — |uOOOOOOd) + |uOOOOdOO)) 

H 24 = ( ) 



Subspace No. 25 with N = 2, s = 0, m s = 0, T 9 _i 



|^25-l) = -7= (|00000002) + |00000020) — |00000200) — |00002000) 

2y 2 

+100020000) + |00200000) — |02000000) — |20000000)) 

1^25-2) = — J (|000000du) — |000000nd) — |0000du00) + |0000ud00) + |000d000u) — |000u000d) 
-f |00d000u0) + |00dn0000) — |00u000d0) — |00ud0000) — |0d000u00) + |0u000d00) 

— |d000n000) — |dn000000) + |u000d000) + |ud000000)) 

1^25-3) = -7= (|000d00n0) — |000n00d0) + |00d0000u) — |00u0000d) 

2v 2 

— |OdOOuOOO) + |OnOOdOOO) — |dOOOOuOO) + |uOOOOdOO)) 

U -2y/2t 

H 25 = ( -2y/2t 4W-^- 2y/2t 
2y/2t 
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Subspace No. 26 with N = 2, s = 0, m s = 0, T 9 _i 

|^26-l) = -7= (|00000002) — |00000020) — |00000200) + |00002000) 

2v 2 

+100020000) — |00200000) — |02000000) + |20000000)) 

1^26-2) = — ^ (|00000dn0) — |00000ud0) — |0000d00u) + |0000u00d) — |000d000u) + |000u000d) 
-f |00d000u0) — |00n000d0) + |OdOOOnOO) + |OduOOOOO) — |OuOOOdOO) — |0nd00000) 

— IdOOOnOOO) — |dOOuOOOO) + |uOOOdOOO) + |uOOdOOOO)) 

1^26-3) = — - 7= (|OOOduOOO) — |OOOudOOO) — |OOdOOuOO) + |OOuOOdOO) 

2y 2 

— |OdOOOOuO) + |OuOOOOdO) + |dOOOOOOu) — |uOOOOOOd)) 

/ U 2V2t 
H 26 = 2y/2t AW 2y/2t 
\ 2y/2t 

Subspace No. 27 with N = 2, s = 0, m s = 0, r 9 _ 3 

1^27-1) = -7= (|00000002) + |00000020) + |00000200) + |00002000) 

2y2 

— |00020000) — |00200000) — |02000000) — |20000000>) 

1^27-2) = — J (|000000du) — |000000ud) + |OOOOOdnO) — |OOOOOudO) + |0000d00u) + |OOOOduOO) 

— |0000n00d) — |OOOOudOO) — |OOduOOOO) + |OOudOOOO) — |OduOOOOO) + |OudOOOOO) 
— IdOOnOOOO) — |duOOOOOO) + |uOOdOOOO) + |udOOOOOO)) 

1^27-3) = -7= (|OOOOOdOn) — |OOOOOuOd) + |OOOOdOuO) — |OOOOuOdO) 

2y 2 

— |OdOuOOOO) + |OnOdOOOO) — |dOuOOOOO) + |uOdOOOOO)) 

/ U -2y/2t 

H 27 = -2y/2t AW-^ 2y/2t 
\ 2V2t 

Subspace No. 28 with N = 2, s = 0, m s = 0, Tio-i 

I ^28-1 ) = — J (|000000dn) — |000000nd) — |OOOOOduO) + |OOOOOudO) — |0000d00u) + |OOOOduOO) 
-f |0000u00d) — |OOOOudOO) — |OOduOOOO) + |OOudOOOO) + |OduOOOOO) — |OudOOOOO) 
-f |dOOuOOOO) — |duOOOOOO) — |uOOdOOOO) + |udOOOOOO)) 

H 28 = ( AW - f ) 
Subspace No. 29 with N = 2, s = 0, m s = 0, Tio-2 

|^29-l) = — J (|000000du) — |000000ud) — |OOOOduOO) + |OOOOudOO) — |000d000u) + |000u000d) 

— |OOdOOOnO) + |OOdnOOOO) + |OOuOOOdO) — |OOudOOOO) + |OdOOOuOO) — |OuOOOdOO) 
-f |dOOOuOOO) — |dnOOOOOO) — |uOOOdOOO) + |udOOOOOO)) 

H 29 = ( AW - f ) 
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Schumann and Zwicker: Extended Hubbard model on a cubic clusters 



Subspace No. 30 with N = 2, s = 0, m s = 0, Tio-3 



| ^30-1 ) = — ^ (|00000du0) — |00000ud0) — |0000d00u) + |0000u00d) + |000d000u) — |000u000d) 
— |00d000n0) + |00n000d0) — |0d000n00) + |0du00000) + |0u000d00) — |0nd00000) 
-f |d000u000) — |d00n0000) — |u000d000) + |u00d0000)) 

H 30 = ( W - f ) 



Subspace No. 31 with N = 2, s = 0, m s = 0, T 7 

|^31_l) = — ^= (|00000002) — |00000020) + |00000200) — |00002000) 

2v 2 

— |00020000) + |00200000) — |02000000) + |20000000)) 

H 3 i = ( U ) 
Subspace No. 32 with N = 2, s = 2, m s = 0, T 4 -i 



| ^32-1 ) = ^ (|000000dn) + |000000ud) — |0000du00) — |0000ud00) — |000d000u) — |000u000d) 

-f |00d000u0) — |00dn0000) + |00u000d0) — |00ud0000) + |0d000u00) + |0u000d00) 

— |d000n000) + |dn000000) — |u000d000) + |ud000000)) 

1^32-2) = J (|00000d0n) + |00000n0d) — |0000d0u0) — |0000u0d0) — |000du000) — |000ud000) 

-f |00d00u00) + |00n00d00) + |0d0000u0) — |0d0u0000) + |0u0000d0) — |0u0d0000) 

— |d000000n) + |d0n00000) — |u000000d) + |u0d00000)) 

Hs2 - I 2t 



Subspace No. 33 with N = 2, s = 2, m s = 0, T 4 _ 2 



|^33-l) = — J (|00000d0n) + |00000n0d) + |0000d0u0) + |0000u0d0) — |000d00u0) — |000u00d0) 

— |00d0000n) — |00n0000d) + |0d00u000) — |0d0u0000) + |0u00d000) — |0u0d0000) 
-f |d0000u00) — |d0n00000) + |u0000d00) — |u0d00000)) 

1^33-2) = — J (|00000du0) + |00000nd0) + |0000d00u) + |0000u00d) — |000d000u) — |000u000d) 

— |00d000n0) — |OOnOOOdO) + |OdOOOuOO) — |OduOOOOO) + |OuOOOdOO) — |OudOOOOO) 
-f |dOOOuOOO) — |dOOnOOOO) + |uOOOdOOO) — |nOOdOOOO)) 

2t 



H 



33 



2t ±(J + 8W) 
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Subspace No. 34 with N = 2, s = 2, m s = 0, T 4 _ 3 



|^34-l) = ^ (|000000dn) + |000000ud) + |00000du0) + |00000ud0) — |0000d00u) + |0000du00) 

— |0000n00d) + |OOOOndOO) + |OOduOOOO) + |OOudOOOO) + |OduOOOOO) + |OndOOOOO> 

— IdOOnOOOO) + |dnOOOOOO) — |uOOdOOOO) + |ndOOOOOO)) 

1^34-2) = — J (|OOOdOOnO) — |OOOdnOOO) + |OOOuOOdO) — |OOOudOOO) — |00d0000u) + |OOdOOuOO) 

— |00n0000d) + |OOnOOdOO) — |OdOOOOuO) + |OdOOuOOO) — |OuOOOOdO) + |OuOOdOOO) 
-f |dOOOOOOu) — |dOOOOnOO) + |uOOOOOOd) — |uOOOOdOO)) 

±(J + 8W) 2t 
2t 



H 



34 



Subspace No. 35 with N = 2, s = 2, m s = 0, T 5 _i 



|^35-l) = J (|000000dn) + |000000nd) — |OOOOOduO) — |OOOOOudO) + |0000d00u) + |OOOOduOO) 

-f |0000u00d) + |OOOOndOO) + |OOduOOOO) + |OOudOOOO) — |OduOOOOO) — |OudOOOOO) 
-f |dOOuOOOO) + |dnOOOOOO) + |uOOdOOOO) + |udOOOOOO)) 

1^35-2) = — J (|OOOdOOuO) + |OOOdnOOO) + |OOOuOOdO) + |OOOudOOO) — |00d0000u) — |OOdOOuOO) 

— |00n0000d) — |OOnOOdOO) + |OdOOOOuO) + |OdOOuOOO) + |OnOOOOdO) + |OuOOdOOO) 

— |dOOOOOOn) — |dOOOOnOO) — |uOOOOOOd) — |uOOOOdOO)) 

_ / ±(J + 8W) 2t \ 
Hs5 - V 2t ) 



Subspace No. 36 with N = 2, s = 2, m s = 0, r 5 _ 2 



|^36-l) = J (|000000du) + |000000ud) — |OOOOduOO) — |OOOOndOO) + |000d000u) + |000u000d) 

— |OOdOOOnO) — |OOdnOOOO) — |OOuOOOdO) — |OOudOOOO) — |OdOOOuOO) — |OuOOOdOO) 
-f |dOOOuOOO) + |dnOOOOOO) + |uOOOdOOO) + |udOOOOOO)) 

1^36-2) = J (|OOOOOdOn) + |OOOOOnOd) — |OOOOdOuO) — |OOOOuOdO) + |OOOduOOO) + |OOOndOOO) 

— |OOdOOnOO) — |OOnOOdOO) — |OdOOOOuO) — |OdOuOOOO) — |OuOOOOdO) — |OuOdOOOO) 
-f |dOOOOOOu) + |dOnOOOOO) + |uOOOOOOd) + |uOdOOOOO)) 



1 {J+8 W) 2t 
Hm ~ 1 2t 
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Schumann and Zwicker: Extended Hubbard model on a cubic clusters 



Subspace No. 37 with N = 2, s = 2, m s = 0, r 5 _ 3 

|^37-l) = — ^ (|00000d0u) + |00000u0d) + |0000d0u0) + |0000u0d0) + |000d00u0) + |000u00d0) 
-f |00d0000u) + |00u0000d) — |0d00u000) — |0d0u0000) — |0u00d000) — |0u0d0000) 

— |d0000n00) — |d0u00000) — |u0000d00) — |u0d00000)) 

1^37-2) = — ^ (|00000du0) + |OOOOOudO) + |0000d00u) + |0000u00d) + |000d000u) + |000u000d) 
-f |OOdOOOuO) + |OOuOOOdO) — |OdOOOuOO) — |OduOOOOO) — |OuOOOdOO) — |OudOOOOO) 

— IdOOOnOOO) — |dOOuOOOO) — |uOOOdOOO) — |uOOdOOOO)) 

it ( 2t 



2t ±(J + 8W) 
Subspace No. 38 with N = 2, s = 2, m s = 0, T 7 

I ^38-1 ) = -7= (|000000du) + |000000ud) — |OOOOOduO) — |OOOOOudO) + |0000d00u) + |OOOOduOO) 

2y6 

-f |0000u00d) + |OOOOudOO) + |000d000u) + |000u000d) — |OOdOOOuO) — |OOduOOOO) 

— |OOnOOOdO) — |OOudOOOO) + |OdOOOuOO) + |OduOOOOO) + |OuOOOdOO) + |OudOOOOO) 

— |dOOOuOOO) — |dOOuOOOO) — |duOOOOOO) — |uOOOdOOO) — |uOOdOOOO) — |udOOOOOO>) 

|^r 38 _ 2 ) = 3— (|OOOdOuOO) + |OOOuOdOO) — |OOdOuOOO) — |OOuOdOOO) 

2y2 

-f |OdOOOOOu) + |OuOOOOOd) — |dOOOOOuO) — |uOOOOOdO)) 

H 38 = ( + 



Subspace No. 39 with N = 2, s = 2, m s = 0, T 8 _i 

|^39_l) = ^= (|000000du) + |000000ud) — |OOOOOduO) — |OOOOOudO) + |0000d00u) + |OOOOduOO) 

4y3 

-f |0000u00d) + |OOOOudOO) — |OOduOOOO) — |OOudOOOO) + |OduOOOOO) + |OudOOOOO) 
— |dOOuOOOO) — |dnOOOOOO) — |uOOdOOOO) — |udOOOOOO)) 

+ (|000d000u) + |000n000d) — |OOdOOOuO) — |OOnOOOdO) 

-f |OdOOOuOO) + |OnOOOdOO) — |dOOOuOOO) — |uOOOdOOO)) 

H 39 = ( W + «w) ) 

Subspace No. 40 with N = 2, s = 2, m s = 0, r 8 _ 2 

I ^40 — 1 ) = — ^ (|000000dn) + |000000ud) + |OOOOOduO) + |OOOOOudO) — |0000d00u) + |OOOOduOO) 
— |0000n00d) + |OOOOndOO) — |OOduOOOO) — |OOudOOOO) — |OduOOOOO) — |OudOOOOO) 
-f |dOOuOOOO) — |dnOOOOOO) + |nOOdOOOO) — |ndOOOOOO)) 
H 40 = ( ^(J + 8W) ) 
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Subspace No. 41 with N = 2, s = 2, m s = 0, T 9 _i 

| ^41-1 ) = J (|00000d0u) + |OOOOOuOd) + |OOOOdOuO) + |OOOOuOdO) — |OOOduOOO) — |OOOudOOO) 

— |OOdOOuOO) — |OOuOOdOO) + |OdOOOOuO) + |OdOuOOOO) + |OuOOOOdO) + |OuOdOOOO) 
-f |dOOOOOOn) + |dOuOOOOO) + |uOOOOOOd) + |uOdOOOOO)) 

1^41-2) = — - 7= (|OOOOOduO) + |OOOOOudO) + |0000d00u) + |0000u00d) 

2y 2 

-f |OduOOOOO) + |OndOOOOO) + |dOOuOOOO) + |uOOdOOOO)) 
1^41-3) = -7= (|OOOdOuOO) + |OOOnOdOO) + |OOdOuOOO) + |OOuOdOOO) 

— |OdOOOOOn) — |OnOOOOOd) — |dOOOOOuO) — |uOOOOOdO)) 

2>/2t 2^2t 

H41 = ( 2y/2t \(J + SW) 
2y/2t 

Subspace No. 42 with N = 2, s = 2, ra s = 0, T 9 _i 

|^42-l) = -7= (|000000dn) + |000000nd) — |OOOOduOO) — |OOOOudOO) 

2y 2 

-f |OOduOOOO) + |OOndOOOO) — |duOOOOOO) — |udOOOOOO)) 

|^42-2) = — J (|OOOOOdOu) + |OOOOOnOd) — |OOOOdOuO) — |OOOOuOdO) — |OOOdOOuO) — |OOOuOOdO) 
-f |00d0000u) + |00n0000d) + |OdOOuOOO) + |OdOuOOOO) + |OuOOdOOO) + |OuOdOOOO) 

— IdOOOOnOO) — |dOnOOOOO) — |uOOOOdOO) — |uOdOOOOO)) 

1^42-3) = — 7= (|OOOdOuOO) + |OOOnOdOO) — |OOdOuOOO) — |OOuOdOOO) 

2y2 

— |OdOOOOOn) — |OnOOOOOd) + |dOOOOOuO) + |nOOOOOdO)) 

/ ±(J + 8W) 2V2t 
H 42 = 2V2t 2V2t 

\ 2y/2t 

Subspace No. 43 with N = 2, s = 2, m s = 0, r 9 _ 3 

|^43-l) = — - 7= (|000d000n) + |000n000d) + |OOdOOOuO) + |OOuOOOdO) 
2y 2 

-f |OdOOOuOO) + |OnOOOdOO) + |dOOOuOOO) + |uOOOdOOO)) 

1^43-2) = J (|OOOdOOnO) + |OOOdnOOO) + |OOOuOOdO) + |OOOudOOO) + |00d0000u) + |OOdOOuOO) 

-f |00u0000d) + |OOnOOdOO) + |OdOOOOuO) + |OdOOuOOO) + |OuOOOOdO) + |OnOOdOOO) 
-f |dOOOOOOu) + IdOOOOnOO) + |uOOOOOOd) + |uOOOOdOO)) 

1^43-3) = — - 7= (|OOOdOnOO) + |OOOnOdOO) + |OOdOuOOO) + |OOuOdOOO) 

2v 2 

-f |OdOOOOOu) + |OnOOOOOd) + |dOOOOOuO) + |uOOOOOdO)) 

±(J + 8W) 2y/2t 
H43 = ( 2y/2t 2y/2t 

2y/2t 
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Schumann and Zwicker: Extended Hubbard model on a cubic clusters 



Subspace No. 44 with N = 2, s = 2, m s = 0, Tiq-i 

| ^44-1 ) = ^ (|000d00n0) — |000du000) + |000u00d0) — |000ud000) + |00d0000u) — |00d00u00) 

-f |00u0000d) — |00u00d00) — |0d0000u0) + |0d00u000) — |0u0000d0) + |0u00d000) 

— |d000000n) + |d0000n00) — |uOOOOOOd) + |uOOOOdOO)) 
H44 = ( ) 

Subspace No. 45 with N = 2, s = 2, m s = 0, Tio-2 

|^45-l) = J (|00000d0n) + |OOOOOnOd) + |OOOOdOuO) + |OOOOnOdO) + |OOOduOOO) + |OOOudOOO) 

-f |OOdOOuOO) + |OOnOOdOO) — |OdOOOOuO) + |OdOuOOOO) — |OuOOOOdO) + |OuOdOOOO) 

— IdOOOOOOn) + |dOnOOOOO) — |uOOOOOOd) + |uOdOOOOO)) 
H45 = ( ) 

Subspace No. 46 with N = 2, s = 2, m s = 0, Tio-3 

|^46-l) = — J (|OOOOOdOu) + |OOOOOuOd) — |OOOOdOuO) — |OOOOuOdO) + |OOOdOOuO) + |OOOuOOdO) 

— |00d0000n) — |00n0000d) — |OdOOuOOO) + |OdOuOOOO) — |OuOOdOOO) + |OuOdOOOO) 
-f |dOOOOuOO) — |dOnOOOOO) + |uOOOOdOO) — |uOdOOOOO)) 

H 46 = ( ) 
Subspace No. 47 with N = 2, s = 2, m s = 1, T4-1 

|^47-l) = — " 7= (|000000nn) — |OOOOuuOO) — |000n000n) + |OOnOOOnO) 

2y 2 

— |OOnnOOOO) + |OuOOOuOO) — |uOOOuOOO) + |uuOOOOOO)) 
1^47-2) = — " 7= (|OOOOOuOu) — |OOOOuOuO) — |OOOuuOOO) + |OOuOOuOO) 

2y 2 

-f |OuOOOOuO) — |OuOuOOOO) — |uOOOOOOu) + |uOuOOOOO)) 

Subspace No. 48 with N = 2, s = 2, m s = 1, r 4 - 2 

|^48-l) = -7= (|OOOOOnOn) + |OOOOuOuO) — |OOOuOOuO) — |00u0000u) 

2y 2 

-f |OuOOuOOO) — |OuOuOOOO) + |uOOOOuOO) — |uOuOOOOO)) 

|^r 48 _ 2 ) = 3— (|OOOOOuuO) + |0000u00u) — |000u000u) — |OOuOOOuO) 

2y2 



H 



48 - + 1 



-f |OuOOOuOO) — |OnnOOOOO) + |nOOOnOOO) — |nOOnOOOO)) 

2t 



2t MJ + 8W) 
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Subspace No. 49 with TV = 2. s = 2. m s = 1. T 4 _ 3 



| ^49-1 ) = — - 7= (|000000uu) + |OOOOOuuO) — |OOOOuOOu) + |OOOOuuOO) 
2y 2 



+ |OOnnOOOO) + |OnnOOOOO) — |nOOnOOOO) + |uuOOOOOO>) 
|^49-2) = -7= (|OOOuOOuO) — |OOOuuOOO) — |OOuOOOOu) + |OOuOOuOO) 

2y2 

— |OnOOOOnO) + |OuOOuOOO) + |uOOOOOOu) — |uOOOOuOO)) 

H49 " 1 2t 



Subspace No. 50 with N = 2, s = 2, m s = 1, T 5 _i 



|^50-l) = — - 7= (|000000nn) — |00000uu0) + |0000u00u) + |0000uu00) 

2y 2 

-f |00uu0000) — |0uu00000) + |u00u0000) + |uu000000)) 
1^50-2) = -7= (|000u00u0) + |OOOuuOOO) — |OOuOOOOu) — |OOuOOuOO) 

2y 2 

+ |OnOOOOnO) + |OnOOnOOO) — |uOOOOOOu) — |uOOOOuOO)) 

/ ±(J + 8W) 2t 
Hso = 1 2t 



Subspace No. 51 with N = 2, s = 2, m s = 1, T5_2 



|^51-l) = — - 7= (|000000nn) — |OOOOnnOO) + |OOOuOOOu) — |OOuOOOuO) 

2y 2 



-|OOuuOOOO) — |OuOOOuOO) + |uOOOuOOO) + |uuOOOOOO>) 
(|OOOOOuOu) — |OOOOuOuO) + |OOOuuOOO) — |OOuOOuOO) 
-|OuOOOOnO) — |OnOnOOOO) + |uOOOOOOu) + |nOnOOOOO)) 



01 " 1 2t 



Subspace No. 52 with N = 2, s = 2, m s = 1, T 5 _ 3 



|^52-l) = -7= (|OOOOOnOn) + |OOOOuOuO) + |OOOuOOuO) + |OOuOOOOu) 

2y 2 

— |OnOOnOOO) — |OuOuOOOO) — |uOOOOuOO) — |uOuOOOOO)) 
1^52-2) = -7= (|OOOOOuuO) + |OOOOuOOu) + |OOOuOOOu) + |OOuOOOuO) 

2y2 

— |OnOOOnOO) — |OuuOOOOO) — |uOOOuOOO) — |uOOuOOOO)) 

_ / 2t 

t±52 - o+l 



It UJ + 8W) 
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Schumann and Zwicker: Extended Hubbard model on a cubic clusters 



Subspace No. 53 with N = 2, s = 2, m s = 1, T 7 



| ^53-1 ) = "7= (|000000nn) — |00000uu0) + |0000u00u) + |0000uu00) + |000u000u) — |00u000u0) 

2y 3 

— |00nn0000) + |0n000n00) + |0uu00000) — |u000u000) — |u00u0000) — |nn000000)) 
1^53-2) = — ^ (|000u0u00) — |00u0u000) + |0u00000u) — |u00000u0)) 

53 ~ 1 



Subspace No. 54 with AT = 2, s = 2, m s = 1, T 8 _i 



|^54-l) = -7= (|000000nn) — |00000uu0) + |0000u00u) + |0000uu00) 

2y6 

— |00nn0000) + |0uu00000) — |u00u0000) — |uu000000)) 
1 / 

+ — (|000u000u) — |00u000u0) + |0u000u00) — |n000n000)) 

V6 

H 54 = ( ^(J + 8W) ) 
Subspace No. 55 with N = 2, s = 2, ra s = 1, r 8 _ 2 



|^55-l) = -7= (|000000nn) + |00000uu0) — |0000u00u) + |0000uu00) 

2y 2 

— |00nn0000) — |OuuOOOOO) + |nOOnOOOO) — |uuOOOOOO)) 

H 55 = ( ±(J + 8W) ) 



Subspace No. 56 with N = 2, s = 2, ra s = 1, r 9 _i 



I ^56-1 ) = — - 7= (|00000u0u) + |OOOOuOuO) — |OOOuuOOO) — |OOuOOuOO) 
2 V 2 

-f |OuOOOOuO) + |OuOuOOOO) + |uOOOOOOu) + |uOuOOOOO)) 
1^56-2) = ^ (|OOOOOuuO) + |OOOOuOOu) + |OuuOOOOO) + |uOOuOOOO)) 

|^ 56 _ 3 ) = — i (|OOOuOuOO) + |OOuOuOOO) — |OuOOOOOu) — |uOOOOOuO)) 
2y/2t 2y/2t 

H 56 = ( 2V2t 1(J + 8W) 
2y/2t ~ 
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Subspace No. 57 with N = 2, s = 2, ra s = 1, T 9 _i 

|^ 57 _l) = — i (|000000uu) — |0000nn00) + |00uu0000) — |uu000000)) 

1^57-2) = -7= (|00000n0n) — |OOOOnOnO) — |OOOnOOnO) + |OOnOOOOn) 

-f |OuOOuOOO) + |OuOuOOOO) — |nOOOOnOO) — |nOnOOOOO)) 
|^ 57 _ 3 ) = - (|OOOnOnOO) — |OOuOuOOO) — |OuOOOOOu) + |uOOOOOuO)) 

/ ±(J + 8W) 2y/2t 
H 57 = 2>/2t 2y/2t 

\ 2V2t 

Subspace No. 58 with N = 2, s = 2, m 8 = 1, r 9 _ 3 

|^r 58 _ 1 ) = 1 (|OOOuOOOu) + |OOuOOOuO) + |OuOOOuOO) + |uOOOuOOO)) 
1^58-2) = — 7= (|OOOuOOuO) + |OOOuuOOO) + |OOuOOOOu) + |OOuOOuOO) 

2y 2 

+ |OnOOOOnO) + |OuOOuOOO) + |uOOOOOOu) + |uOOOOuOO)) 
|^ 58 _ 3 ) = i (|OOOuOuOO) + |OOuOuOOO) + |OuOOOOOu) + |uOOOOOuO)) 

\(J + SW) 2V2t 
H 58 = ( 2y/2t 2y/2t 

2V2t 

Subspace No. 59 with N = 2, s = 2, ra s = 1, Tiq-i 

|^59_l) = —^7= (|OOOuOOuO) — |OOOuuOOO) + |OOuOOOOu) — |OOuOOuOO) 
2y 2 

— |OnOOOOnO) + |OnOOnOOO) — |nOOOOOOn) + |uOOOOuOO)) 
H59 = ( ) 

Subspace No. 60 with N = 2, s = 2, ra s = 1, Tio-2 

I ^60-1 ) = ^= (|00000n0n) + |0000u0u0) + |000uu000) + |00u00u00) 

— |0n0000n0) + |0u0u0000) — |u000000u) + |u0u00000)) 

H 60 = ( ) 
Subspace No. 61 with N = 2, s = 2, m s = 1, Tio-3 

I ^61-1 ) = -7= (|00000n0n) — |0000u0u0) + |000u00u0) — |00u0000u) 

2y 2 

— |0n00n000) + |0n0n0000) + |u0000u00) — |u0u00000)) 

H 6 i = ( ) 
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Schumann and Zwicker: Extended Hubbard model on a cubic clusters 



B Groundstate of the Hubbard model for U = oo and N e = 6 



In the following we give the groundstate of the pure Hubbard model with infinite on-site interaction 
U and six electrons, i.e. two holes in the half filled cluster. 



\9 G s) = \s = 0,m a = 0,r 1 ,E = -2y/4+V3t) 

= Ci|$i) + C 2 |$ 2 > + C 3 |$ 3 > + C 4 |$ 4 > 

+ C 5 |$ 5 ) + C 6 |$ 6 ) + C 7 |$ 7 > + ^ 8 |* 8 > + Cq |$ 9 > 



(9) 



$l) = (\00ddduuu) — \00dduduu) — \00dudduu) — \00duuudd) + |00ndddnn) + \00uduudd) 



+ 


OOuududd) — 


\00uuuddd) + 


\0dd0duuu) — 


\0dd0uuud) — 


\0ddu0udu) -\- 


\0du0duud) 


+ 


OduOuddu) — 


\0dud0duu) — 


\0dud0uud) -\- 


\0duu0dud) — 


\0ud0duud) — 


\0ud0uddu) 




OuddOudu) -\- 


\0udu0ddu) -\- 


\0udu0udd) -\- 


\0uu0dddu) — 


\0uu0uddd) -\- 


\0uud0dud) 


+ 


dOOduduu) — 


\d00duudu) -\- 


\d00uduud) -\- 


\d00uuddu) — 


\d0dud0uu) — 


\d0duu0du) 




dOuduOud) -\- 


\d0uud0du) -\- 


\dd00uudu) — 


\dd00uuud) — 


\dd0udu0u) -\- 


\ddduOuuO) 


+ 


ddduuuOO) — 


\dduOuduO) -\- 


\ddudu00u) — 


\dduduu00) — 


\dduu00du) -\- 


\dduu00ud) 




dduuduOO) -\- 


\dduuud00) — 


\du00dduu) — 


\du00uudd) — 


\du0ddu0u) -\- 


\du0udd0u) 


+ 


duOuduOd) — 


\dudOduuO) — 


\dudOuudO) -\- 


\dudd00uu) — 


\duddu00u) -\- 


\duuOdudO) 


+ 


duudOduO) — 


\duudOudO) -\- 


\duudd00u) — 


\duudu00d) -\- 


\duuu00dd) -\- 


\duuuOddO) 




uOOdduud) — 


\u00duddu) -\- 


\u00uddud) — 


\u00ududd) — 


\u0ddu0ud) -\- 


\u0dud0du) 


+ 


uOuddOud) -\- 


\u0udu0dd) -\- 


\ud00dduu) -\- 


\ud00uudd) — 


\ud0dud0u) — 


\ud0duu0d) 


+ 


udOuudOd) — 


\uddOuduO) — 


\uddd00uu) — 


\udddOuuO) -\- 


\udduOduO) — 


\udduOudO) 


+ 


uddudOOu) — 


\udduu00d) -\- 


\uduOdduO) -\- 


\uduOuddO) — 


\uduu00dd) -\- 


\uduud00d) 


+ 


uuOOdddu) — 


\uu00ddud) -\- 


\uu0dud0d) -\- 


\uudOdudO) — 


\uudd00du) -\- 


\uudd00ud) 




uuddduOO) -\- 


\uuddud00) — 


\uudud00d) -\- 


\uududd00) — 


\uuudOddO) — 


\uuuddd00)^) 



| $2) = (\00dduudu) — \00dduuud) + \00uudddu) — \00uuddud) — \0dd0uduu) + \0dd0uudu) 

-\-\0ddu0duu) — \0duu0ddu) -\- \0udd0uud) — \0uu0ddud) -\- \0uu0dudd) — \0uud0udd) 

— \ d00dduuu) -\- \d00duuud) + \d0udd0uu) — \d0uud0ud) -\- \dd00duuu) — \dd00uduu) 
-\- 1 ddOuudOu) -f- \dddu00uu) — \ddduu00u) -\- \dduOduuO) — \ddud00uu) — \ddudOuuO) 
-\- 1 duOduuOd) -\- \duddOuuO) -\- \dudduu00) — \duuOdduO) — \duuud00d) -\- \duuudd00) 

— \u00udddu) -\- \u00uuddd) -\- \u0ddu0du) — \u0duu0dd) — \ud0udd0u) -\- \uddOuudO) 
-\-\udddu00u) — \uddduu00) — \uduuOddO) — \uduudd00) -\- \uu00dudd) — \uu00uddd) 

— \uu0ddu0d) — \uudOuddO) -\- \uudu00dd) -\- \uuduOddO) — \uuud00dd) -\- \uuudd00d)^) 
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1^3) = (\00dududu) — \00duduud) + \00ududdu) — \00ududud) + \OdduOuud) + |OdduduOu) 



— 


OdduuuOd) — 


\OduOudud) -\- 


OduOuudd) -\- 


OdududOu) -\- 


OduduuOd) — 


\OduuOudd) 


+ 


OduuduOd) — 


\OduuudOd) — 


OudOdduu) -\- 


OudOdudu) -\- 


OuddOduu) -\- 


\OuddduOu) 


— 


OuddudOu) — 


\OududdOu) — 


OududuOd) — 


OuudOddu) -\- 


OuudddOu) — 


\OuududOd) 


— 


dOOududu) -\- 


\dOOuuudd) -\- 


dOduduuO) -\- 


dOduuudO) -\- 


dOuduOdu) -\- 


\dOududuO) 


— 


dOuduudO) — 


\dOuududO) — 


dOuuuOdd) -\- 


dOuuuddO) -\- 


ddOuOudu) — 


\ddOuOuud) 


+ 


ddOuuuOd) — 


\dduOuOdu) -\- 


dduOuOud) -\- 


dduOuudO) — 


dduuOudO) — 


\dduuuOOd) 


+ 


duOOdudu) — 


\duOOuddu) — 


duOdOduu) -\- 


duOdOudu) -\- 


duOdudOu) — 


\duOuOddu) 




duOuOudd) -\- 


\dudOdOuu) -\- 


dudOuOdu) -\- 


duduOOdu) -\- 


duduOudO) — 


\dududOOu) 


+ 


dududuOO) — 


\duuOdOdu) -\- 


duuOuOdd) — 


duuOuddO) — 


duudOOdu) -\- 


\duududOO) 




uOOddduu) -\- 


\uOOdudud) -\- 


uOdddOuu) — 


uOddduuO) -\- 


uOdduduO) — 


\uOdudOud) 


+ 


uOdudduO) — 


\uOdududO) — 


uOuddduO) — 


uOududdO) -\- 


udOOduud) — 


\udOOudud) 


+ 


udOdOduu) -\- 


\udOdOuud) — 


udOuOdud) -\- 


udOuOudd) — 


udOuduOd) — 


\uddOdOuu) 


+ 


uddOduuO) -\- 


\uddOuOud) -\- 


udduOOud) — 


udduduOO) — 


uduOdOud) — 


\uduOuOdd) 




ududOOud) — 


\ududOduO) -\- 


ududuOOd) — 


udududOO) -\- 


uuOdOddu) — 


\uuOdOdud) 




uuOdddOu) — 


\uudOdOdu) -\- 


uudOdOud) — 


uudOdduO) -\- 


uuddOduO) -\- 


\uudddOOu)^) 



I $4) = (\00duuddu) — \00udduud) + \OdddOuuu) — \OduOdduu) + \OudOuudd) — \OuuuOddd) 

— \dOOudduu) -\- \dOdduOuu) -\- \ddOduuOu) -\- \dddOuuuO) — \dduuOduO) — \dduudOOu) 
-\-\duOOduud) — \duudOOud) -\- \duudduOO) -\- \uOOduudd) — \uOuudOdd) — \udOOuddu) 
-\-\udduOOdu) — \udduudOO) — \uuOuddOd) -\- \uuddOudO) -\- \uudduOOd) — \uuuOdddO)^) 

1^5) = (\OdOdduuu) — \OdOduudu) — \OdduuduO) — \OduddOuu) — \OdudduuO) + \OduudOdu) 

-\-\OuOuddud) — \OuOuuddd) — \OudduOud) -\- \OuduuOdd) -\- \OuduuddO) -\- \OuuddudO) 
-\-\dOdOuduu) — \dOdOuuud) — \dOduOduu) — \dOduudOu) — \dOudduOu) -\- \dOuuOdud) 

— \ ddOuuduO) -\- \ddduuOuO) — \dduOduOu) -\- \ddudOuOu) — \duOduOud) -\- \duOudOud) 
-\-\duOudduO) — \dudOOuud) — \dudOuuOd) — \dudduOuO) -\- \duuOOdud) -\- \duuuOdOd) 

— \uOddOudu) -\- \uOduudOd) -\- \uOuOdddu) — \uOuOdudd) -\- \uOudOudd) -\- \uOudduOd) 

— \udOduOdu) — \udOduudO) -\- \udOudOdu) — \uddOOudu) — \udddOuOu) -\- \uduOOddu) 
-\-\uduOddOu) -\- \uduudOdO) -\- \uuOddudO) -\- \uudOudOd) — \uuduOdOd) — \uuuddOdO)^) 

\$>6) = (\OdOudduu) — \OdOuuddu) — \OddduOuu) — \OddduuuO) + \OdduuOdu) — \OduudduO) 

-\-\OuOdduud) — \OuOduudd) -\- \OudduudO) — \OuuddOud) -\- \OuuudOdd) -\- \OuuudddO) 

— \ dOddOuuu) — \dOdduuOu) + \dOuOdduu) — \dOuOduud) -\- \dOudOuud) — \dOuuddOu) 

— \ ddOduOuu) — \ddOduuuO) -\- \ddOudOuu) — \dddOOuuu) — \dddOuuOu) -\- | dduOOduu) 
-\-\dduuOdOu) -\- \dduudOuO) -\- \duOdduuO) — \duuOduOd) -\- \duudOuOd) — \duuddOuO) 
-\-\uOdOuddu) — \uOdOuudd) -\- \uOdduuOd) — \uOduOddu) -\- \uOuuOddd) -\- \uOuuddOd) 

— \udOuuddO) -\- \uddOudOu) — \udduOdOu) -\- \udduuOdO) — \uuOduOdd) -\- \uuOudOdd) 
-\-\uuOudddO) — \uudOOudd) — \uuddOuOd) — \uudduOdO) -\- \uuuOOddd) -\- \uuuOddOd)^) 
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Schumann and Zwicker: Extended Hubbard model on a cubic clusters 



|^>7) = (\0d0uduud) — \0d0uuudd) — \0dduduu0) — \0dudu0ud) — \0dududu0) + \0duuu0dd) 



-\-\0u0ddduu) — \0u0duddu) — \0uddd0uu) -\- \0udud0du) -\- \0ududud0) -\- \0uududd0) 

— \ d0du0udu) — \d0dudu0u) -\- \d0u0uddu) — \d0u0uudd) — \d0udud0u) -\- \d0uu0udd) 

— \ dd0uuud0) — \ddu0uu0d) -\- \dduu0u0d) -\- \dduuu0d0) — \du0du0du) -\- IduOudOdu) 
H~|(iu0ndnd0) — \dud00udu) — \dud0du0u) -\- IduuOOddu) + IduudOdOu) — |duuduO<iO) 
-\-\u0d0dduu) — \u0d0duud) — \u0dd0duu) -\- |n0dwdn0d) + |n0nd0dnd) + \u0udud0d) 

— \ud0du0ud) — \udOduduO) -\- IndOndOnd) — |ndd00nnd) — \uddu0u0d) -\- |ud<iudO«0) 
-h|ndn00dnd) + |ndn0wd0d) + |nnOdddnO) + |u«dO<i<iO«) — \uudd0d0u) — \uudddOuO)^) 



\<&s) — (\0ddduu0u) — \0dduud0u) + \0duudd0u) — \0udduu0d) + |0nnddn0d) — |0uimdd0<i) 

H-|(iOddnnnO) — |dOudduuO) + \dOuudduO) -\- |ddO<iOnnn) — IddOuOduu) + | dddOuOuu) 
— \ddu0d0uu) — \du0d0uud) -\- |duu0d0u<i) — |uOd<iuu<iO) + |nOdnnddO) — \uOuudddO) 



|3>g) = (|0duddu0u) — \OuduudOd) + |dOduuduO) — |duOuOdud) + |dudOuOud) — |nOnddndO) 
-|-|ndOdOndn) — \uduOdOdu)^j 



-\-\udOuOddu) — \uddOuOdu) -\- |uuOdOudd) — |nitOitOddd) + \uudOuOdd) — \uuuOdOdd)^j 



Ci 



J\ (4+^(3 + 2^/3) 



C 2 



180 + 84 V3 
1 





60 + 28V3 



d = 



V /|(4 + V3)(3 + 2V3) 



90 + 42\/3 



C 5 = - 



1 



12V2 




c 8 = 



51 + 2973 




Greek symbols — w-greek . sty 



ON 

o 
o 

a 

On 



-I— > 

e 

O 

o 



> 

On 
00 

in 

m 

o 

On 
O 





a \alpha 




6> \theta 








T 


\tau 




/? \beta 




\vartheta 


7T 


\pi 


V 


\upsilon 




7 \gamma 




/ \ i n1~ a 


TJJ 


\vairpi 


6 


\Dhi 




5 \delta 




/ V \ 11 CJ. ky ky \JL 


n 
H 


\rho 


CO 


\ Vr)TT)h i 




e \epsilon 




A \ lambda 


Q 


\varrho 


X 


\chi 




e \varepsilon 


/x \mu 


a 


\sigma 




Ypsi 




( \zeta 




v \nu 




\var sigma 


00 


\omega 




rj \eta 




£ \xi 












r \itGamma 




yl \itLambda 


E 


\itSigma 




\itPsi 




A \itDelta 




S \itXi 


r 


\itUpsilon 


Q 


\itOmega 




\itTheta 




77 \itPi 




\itPhi 












Table 1: Slanted greek letters 






a 


\upalpha 


e 


\uptheta 





\upo 


X 


\uptau 


p 


\upbeta 




\upvartheta 


71 


\uppi 


V 


\upupsilon 


If 


\upgamma 


i 


\upiota 


05 


\upvarpi 


Y 


\nnnh i 

\ Ulky k^J.±_L 


5 


\updelta 


K 


\upkappa 


P 


\uprho 


(D 

T 


\ uDvarDh i 


c 
c 


\upepsilon 


X 


\uplambda 


Q 


\varrho 


X 


\ n V» n 

\ upciii 


e 


\varepsilon 




\upmu 


G 


\upsigma 


¥ 


\uppsi 




\upzeta 


V 


\upnu 




\upvar sigma 


CO 


\upomega 




\upeta 




\upxi 










r 


\Gamma 


A 


\Lambda 


E 


\Sigma 




\Psi 


A 


\Delta 




\Xi 


T 


\Upsilon 


n 


\0mega 


e 


\Theta 


n 


\Pi 


$ 


\Phi 







Table 2: Upright greek letters 



1 



LX 


\ "Km 4 \ c» H t^Vi o >- 
\ Dill i_ \d.±pilci_r 


C7 


\Ulll\ \ L^IlcX»cL_r 


LP 


\"hm-Tnl- 

\ uiii\t)_r 


T 


\ "Km -T \ +■ a i i \ 


n 
p 


\ Ulll \ \ Ut3 U a/ 




\ V\tyi -4^ \ t"I~Vi ofal- 
\ Ulll \ \ V CLJ. UlltJUCLj 


IX 


\U111\ \p±j 


1 7 

U 


\ Kityi-TX iitvc; "i 1 nn )■ 
\ Ulll \ \ LLUoJ-XUIlj 


7 


\ Km -T \ rr^TTTm V 
\ Ulll L \HjCUlllllcLj 




\ \\xc\ -T \ "i ri1"al- 
\ Ulll \ \1UU CLj 


TT7 

UJ 


\ KlTYl -T \ T"JT\ "1 \ 
\ Ulll \ \ V CL± U _L J 




\U111\ \pn±j 




\bm{\delta} 




\bm{\kappa} 


P 


\bm{\rho} 




\bm{\varphi} 


e 


\bm{\epsilon} 


A 


\bm{\ lambda} 


Q 


\bm{\varrho} 


X 


\bm{\chi} 




\ Ulll \ \ V cLX c3po±±Ull f 




\ "Km -T \ m n J- 
\ Ulll \ \111U. f 


U 


\ Km -T \ a "i rrm o V 
\ Ulll \ \o±H > llld._r 


w 


\ Km -T \ t\ c "i V 

\uini \psi j 


c 
s 


\ Km -i \ t q +■ o J- 
\ Ulll \ \Z.c3 UcLj 


Is 


\ "Km -T \ m l J- 
\ Ulll \ \I1U. J 


s 


\ Km -I \ trove "i rrm c\ V 
\ Ulll \ \ V dX o -L^lHa. f 


1 . 1 


X "Km -I X ATTlQfTSl >- 

\ Ulll l XUlllcfc^cLj 


?7 


\Dnii. \eta_r 




\Dmi\xi_r 










r 


\bm{\itGamma} 




\bm{\itLambda} 


£ 


\bm{\itSigma} 


& 


\bm{\itPsi} 




\bm{\itDelta} 




\bm{\itXi} 


Y 


\bm{\itUpsilon} 


Q 


\bm{\itOmega} 





\bm{\itTheta} 


iT 


\bm{\itPi} 


& 


\bm{\itPhi} 







Table 3: Boldface variants of slanted greek letters 



a 


\pmb{\upalpha} 


e 


\pmb{\uptheta} 





\pmb{\upo} 




\pmb{\uptau} 


P 


\pmb{\upbeta} 




\pmb{\upvartheta} 


% 


\pmb{\uppi} 


1) 


\pmb{\upupsilon} 


Y 


\pmb{\upgamma} 


i 


\pmb{\upiota} 


rs 


\pmb{ \up varp i } 




\pmb{\upphi} 


8 


\pmb{\updelta} 


K 


\pmb{\upkappa} 


P 


\pmb{\uprho} 




\pmb{\upvarphi} 


e 


\pmb{\upepsilon} 


X 


\pmb{\uplambda} 


Q 


\pmb{\varrho} 


X 


\pmb{\upchi} 


£ 


\pmb{\varepsilon} 




\pmb{\upmu} 


G 


\pmb{\upsigma} 




\pmb{\uppsi} 


c 


\pmb{\upzeta} 


V 


\pmb{\upnu} 




\pmb{\upvarsigma} 


CO 


\pmb{\upomega} 


n 


\pmb{\upeta} 




\pmb{\upxi} 










r 


\bm{\Gamma} 


A 


\bm{\Lambda} 


s 


\bm{\Sigma} 




\bm{\Psi} 


A 


\bm{\Delta} 


i— i 


\bm{\Xi} 


T 


\bm{\Upsilon} 




\bm{\0mega} 





\bm{\Theta} 


n 


\bm{\Pi} 




\bm{\Phi} 







Table 4: Boldface variants of upright greek letters 



2 



